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Abstract 

Detailed fluctuation theorem, a microscopic version of the steady state fluctuation theorem, 
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has been proposed by Jarzynski and demonstrated in the case of Hamiltonian systems weakly 
coupled with reservoirs. We show that an identical theorem for phase space compressibility rate 
can be derived for systems evolving under non-Hamiltonian extended system dynamics, without 



certain limiting assumptions made in the original work. Our derivation is based on the non- 
O 

Hamiltonian phase space formulation of statistical mechanics and does not rely on any assumptions 

a . 

of thermodynamic nature. This version of the detailed fluctuation theorem is shown to be generic 
enough to be applicable to several thermostatting schemes. It is shown that in equilibrium, this 
detailed fluctuation theorem boils down to the detailed balance equation and it is further shown to 

a. 

i , reproduce the Jarzynski's work theorem for driven systems. 
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I. INTRODUCTION 

A family of remarkable results in non-equilibrium statistical mechanics, collectively called 
Fluctuation theorems, has been obtained in recent years. These theorems address the issue 
of how macroscopic irreversibility arises from microscopic time reversible dynamics, and in 
a sense quantify the probability of observing Second law violating events. These fluctu- 
ation theorems are valid for systems arbitrarily far from thermodynamic equilibrium and 
have been demonstrated for both deterministic and stochastic evolution. It is this sweep- 
ing generality that make these results quite extraordinary. One version of the fluctuation 
theorem was first discovered by Evans and Searles in the context of molecular dynamics 
simulation J8j of a steady state system subject to a Gaussian isokinetic thermostat. The 
system under consideration was a steady state system subjected to external forcing. Under 
such circumstances, stationarity can be attained only under constant dissipation of energy. 
In a dynamical system, this dissipation manifests itself as a contraction of the phase space 
volume. Considering a set of phase points centered around Tq(pq, go) at time t=0, under the 
dynamics, we have the mapping dFo— ^dl^ with: 



dTc. 



dT, 



where 



exp(- f A(T s )ds) (1) 



A(T.) = ~f t (2) 



then the heat exchange rate with the thermostat, at temperature T, is 

Q(T,t) = -TA(T,t) (3) 

which provides the link between phase space compression and entropy production rate. 
Evans and Searles then define a dissipation function : 



(4) 



fl t (Xo) = In 


r/(r o)" 

L/(n,o)J 


t 

- [A(T s ,s)ds 




where /(r s , s) is the phase space probability density. 


For a dissipative system, 




w> 






(5) 



When the initial probability density /(r ,0) is drawn from an equilibrium distribution and 
the system is driven away from it by external forcing, the statement of the Evans-Searles 
(transient) fluctuation theorem takes the form: 

mt = A) exp(At) (6) 



P(n t = -A) 

In the case where the initial phase space points are drawn from a steady state distribution, 
Evans and Searles derived a rearranged form of the above equation called the steady state 
fluctuation theorem valid for the long time limit, given by the expression: 

Um -In P ^ = A) =A (7) 

t-+oct P(tt = -A) K ' 

Independent of this work, Galavotti and Cohen derived the steady state fluctuation theorem 
using average compression as: 

,. 1, P(-R = A) 

hm-ln—^ -r = A (8) 

These fluctuation theorem are valid in the non-linear, nonequilibrium regime where very 
general results are available. Fluctuation relation have been derived under a variety of 
conditions [[8]] and have been demonstrated experimentally in small systems. 

II. JARZYNSKI'S DETAILED FLUCTUATION THEOREM 



In an interesting development, Jarzynski 
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2] derived a hybrid fluctuation theorem, which 



is also a statement of detailed balance. This detailed fluctuation theorem shows that under 
a nonequilibrium process, the ratio of the probability of observing specific trajectory - anti- 
trajectory pairs goes as an exponential of the entropy produced. For a trajectory that starts 
in a microstate za and evolves into zb in a duration r, resulting in the production of entropy 
AS, its anti-trajectory is the one that starts from z B and evolves into z A causing an entropy 
consumption of AS" . Here (q, p)* stands for reversal of momenta (q, — p) . The detailed 
fluctuation theorem, like all fluctuation theorems, shows that it is more likely that entropy 
is generated rather than consumed in a non-equilibrium process. This theorem is however 
distinct from the other forms of fluctuation theorems in that it makes specific reference to 
the initial and final microstates. 



Jarzynski has derived this result for a Hamiltonian system weakly coupled to a set of 
Hamiltonian reservoirs. The system is manipulated by an external protocol which involves 
making or breaking contact with external forces and heat reservoirs in a specified sequence. 
In this scheme, the application of the external protocol in the reversed order amounts to the 
realization of the time-reversed trajectories. A process n + is defined to be the execution 
of a given protocol for a time interval r and correspondingly, the process II - is the same 
protocol executed in the reverse order. The detailed fluctuation theorem then assumes the 
form 

P '"" A * 1 *' 1 e AS/k B (g) 



P+(z B ,AS\z A ) ^ AS/k£ 



P-(z* A ,-AS\z* B ) 

where the numerator (denominator) in LHS is the joint conditional probability that, under 
the process n + (n~), that the system starts from za {z b ) and ends at zb (z a ), producing 
(consuming) an entropy AS. The derivation of the theorem follows directly from the as- 
sumptions that the reservoir temperatures do not change and their degrees of freedom at 
both initial and final times are Maxwell-Boltzmann distributed. The entropy generated by 
the dynamical evolution is then assumed to be: 

As = -jr^h (io) 

71=1 n 

where AQ n is the the change in the internal energy of the n th heat reservoir. These as- 
sumptions, coupled with the fact that Hamiltonian evolution is time reversible and phase- 
space conserving, lead directly to the detailed fluctuation theorem. Regarding the particular 
choice made for AS, eq lfTOj) . Jarzynski explicitly states that it is valid only when the system 
is weakly coupled with the heat reservoirs possessing infinite number of degrees of freedom. 
Towards the end of his paper, Jarzynski suggests that it would be desirable to perform 
experiments to test the detailed fluctuation theorem. This requires experimental control 
and manipulation of the full microstates of the system which is very difficult to realize ex- 
perimentally. This might explain only very few results are available on the experimental 
verification of the detailed fluctuation theorem in contrast to work fluctuation theorem of 
Jarzynski. In contrast to the experimental scenario, in molecular dynamics one has full 
control and information about the specific microstate and complete dynamics of the system 
and hence one can attempt to verify the detailed fluctuation theorem in such a setting. 



III. DETAILED FLUCTUATION THEOREM IN MOLECULAR DYNAMICS 

In the context of molecular dynamics simulations, thermostatting is realized by appending 
a few extra degrees of freedom to the system, the time evolution of which is of a non- 
Hamiltonian character and hence cannot be understood as a simple perturbation of the 
original uncoupled dynamics. Hence, the question of whether these fluctuation theorems 
can be captured in molecular dynamics simulations has been attracting a lot of attention 
in recent years. For example, the work fluctuation theorem has been rederived for Nose- 
Hoover dynamics without the weak coupling assumption [9|. Similarly, the validity of these 
luctuation theorems for different molecular dynamics ensembles has appeared in literature 
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13]. Within the context of detailed fluctuation theorem, Jarzynski |l|, |2j, has suggested 
that it would be interesting to see whether the DFT manifests itself in the Nose-Hoover 
thermostatted dynamics. Nose-Hoover dynamics belong to a class of thermostatting schemes 
called Extended Phase Space Methods [3J and we realize that in the case of such systems, 
there are no coupling terms in the total Hamiltonian and hence one cannot talk of a weak 
coupling or even delineate the system energy from the heat bath energy, as was done in the 
original derivation. Therefore one cannot use the canonical definition of entropy (eq 10) . 
Further, in the extended system dynamics, temperature enters only as an external parameter 
and hence taking recourse to thermodynamic definition of entropy or internal energy is not 
viable. Given a set of autonomous equations of motion, perhaps the only quantity that can 
be related to entropy is the phase space compressibility. Although the absolute entropy for 
such systems may not be definable, one can identify, under certain assumptions, phase space 
compression rate with the entropy production rate 5K7|. Independent of whether such an 
identification is to be made or not, we demonstrate that a detailed fluctuation theorem can 
be formulated for the phase space compressibility that takes exactly the same form as that 
of the original. 

IV. DERIVATION OF THE DETAILED FLUCTUATION THEOREM FOR THE 
PHASE COMPRESSION RATE 

The present derivation follows essentially the same methodology as that of the original, 
but differs significantly in that no use is made of the Liouville theorem since, in general, 



phase space volume is not conserved under a non-Hamiltonian evolution. The notations and 
symbols used are identical to the ones used in the original, so as to make the similarities 
and distinctions between the two derivations easily observable. 

We start by considering a system ip coupled to a thermostat such that, the extended 
system undergoes a time-reversible, non-Hamiltonian dynamics for which a constant energy 
functional H, hence-forth referred to as the psuedo-Hamiltonian, can be identified. It is 
possible to cast many of the present day thermostats and ergostats into this form, see (4|. 
Let z = (r, p) denote the phase space of the system ip and let the degrees of freedom of the 
thermostat be represented by Y. Let V = (z, Y) denote a point in the phase space of the 
extended system, evolving under the deterministic non-Hamiltonian dynamics: |4| 

2Nd „ 

r . = r B ~ (ii) 

for i = 1, ..., 2Nd, where "H(T) is the generalized energy and B is a 2Ndx2Nd anti-symmetric 
matrix. The equations are known to generate energy preserving dynamics. For instance, 
1 

-1 

ate dimensions, these equations reduce to the usual Hamilton's equations of motion. Any 

other choice of matrix B, so long as it is anti-symmetric, still generates a non-Hamiltonian 
dynamics that conserves H: 

i=l j=l J 

We further note that time reversibility of the underlying dynamics enforces the condition 



when B 



where and 1 indicate the zero and identity matrices of appropri- 



n(T) = H(T*) (13) 

where T* is obtained from T by reversing all the momentum-like variables. The phase 
space contraction rate A(r s ) at time t = s, defined earlier, may be written as 






We now define entropy production in a time interval r under a non-Hamiltonian evolution 

be 



AS = k B A(r)dt (15) 

o 
where k B denotes the Boltzmann constant. For a given state T = (z, Y), we define a time 

reversed state T* to be (z*, Y*) . Let T + (t) denote the state of the extended system at 

time t starting from an initial extended microstate T° and let T_(£) denote state of the 

system at time t but evolving from the momentum reversed final state F T *. With this as 

the background, we are interested in computing the probability that the system ip, at t — 0, 

starts in a particular microstate z(0) = z A and after evolving for a time t = r, reaches 

a state z(r) = z B , generating an entropy AS. Note that the evolution of F = (z,Y) is 

deterministic, eq ffTTl) . We resort to probabilistic description when we are focusing only on 

the evolution of z. For a particular choice of the initial condition z = z Al the probability 

distribution is over the initial conditions of heatbath variables Y 

P + (z B , AS\z A ) = JdYp(Y)5[z B - z T + (z A ,Y)].5[AS - AS(z A ,Y)] (16) 

where z^(za,Y) and AS'(za,Y) denote the final state of the system if) and the net entropy 
produced respectively, after ip has evolved for a time t — r starting at z A with the reservoir 
starting at Y. Using the identity za = z°(z,Y) we can now integrate over the full phase 
space vector V = (z, Y): 

P + (z B , AS\z A ) = J dTp(T)5[z A - z° (r)].<J[z fl - zl(T)].5[AS - AS(T)} (17) 

if the dynamical system eq fill I) possess n c independent conservation laws: 



A k (T) = X k , k = l,..,n e (18) 

then the probability distribution of the mirostates can be written as: [3|] 



p(r) = l[5(A k (T) - X k ) (19) 

fc=i 

In instances where "H is the only conserved quantity of the dynamics, the states T are those 
drawn from the constant H, surface 



p(r°) = 8{H{T°) - E) (20) 

hence, eq (fT7|) may be rewritten as 

P + (z B , AS|z A ) = J dT°5(H(T°) - £)5[z A - z (r°)].5[zB - z T (T°)].S[AS - AS(T )} (21) 
As the dynamics is time reversal invariant, and "H is conserved and we have: 

AS(T T *) = -AS(T°) (22) 

U(T°) = U(T T *) (23) 

z°(r°) = [z T (T T *)]* (24) 
With these expressions at hand, we may now recast the above integral as: 

P + (z B , AS\z A ) = f dT°5(H(T T *) - E)S[z A * - z r (r T *)].<5[z^ - z°(T T *)].5[AS + AS(r T *)] 

(25) 
where we have cast all the variables in the r.h.s of the above integral with the time reversed 
counterparts of the final states, except for the integration volume element dT°. In general, 
as the phase space volume is not conserved in a non-Hamiltonian evolution, dT° is not 
an invariant volume under eq lTTTj) and hence <fT°cannot be directly replaced by dT T * . The 
Jacobian of transformation from the initial to the final time reversed phase space coordinates 
J(T T *; r°) is not unity and we have 

dT T * = J(T T *; T°)dT° = J(T T *; T T )J(T T ; T°)dT° (26) 

As the dynamical system eq (lll|) is assumed to be time reversal invariant, 

j(r r * ; r r ) = i (27) 

therefore, 

j(r T * ; r°) = j(r r ; r°) (28) 



:S 



The time evolution of Jacobian of transformation is 

d_j _ fimm 

The sum on the right is nothing but the negative of the phase space compressibility rate 
A(r 4 ) , hence: 

- fA(r*)dt A~s(r°) 

J(r T ; r°)=e « = e "^^ (30) 

From eqQ22), we know that AS{T T *) = -AS{T°) and we have, 

A"s(r T ') 

dT T * = e~^^dT° (31) 

We may now replace dT in eq(|25]) with e fc s dT T * and we have: 



P+(z B ,AS\z A )= /"e- i ^^rfr T *,5(H(r r *)- J E) ( 5[z A -z (r r *)].5[z B -z T (r")]. ( 5[A 1 S+A5(r r *)], 

(32) 

AS 

We can pull out the factor e fe s out of the integration due to the presence of the delta function 
8[AS + AS(T T % 



P + (z B ,AS\z A ) = e^i f dT T *8{n{T T *)-E)6[z A -z°{T T %S[z B -z T (T T *)].S[AS+AS{T T % 

(33) 

AS 

The right hand side of the above equation is nothing but e fc s times the probability distri- 
bution P(z* B , — AS\z* A ) . Thus we have the final result of the Detailed Fluctuation theorem: 



P(z B , AS\z A ) = e^P(z B , -AS\z* A ), (34) 

In summary, from the ingredients of the above derivation we list all the requirements that 
a dynamical system has to satisfy for the present derivation to go through: 

1 . The phase space of the extended system is of even number of dimensions. 

2. The evolution is deterministic and equations of motion are time reversal invariant. 

3. The phase space compressibility of the system is non-zero. 



4. A constant of motion, usually a "psuedo-Hamiltonian" can be identified for the system 
that is time reversal invariant. 

From the above points, it can be see that no assumptions have been made about the 
specific type of interactions present in the system and there are no restrictions on system 
size and time scale r over which these fluctuations are realized. It should also be appreci- 
ated that this detailed fluctuation theorem is valid arbitrarily far from equilibrium, as no 
assumptions of thermodynamic nature has been made in the derivation. Further, unlike in 
the original derivation, there is no requirement for the system to be driven out of equilib- 
rium through a given protocol. As there are fluctuations in the phase space compressibility 
even at equilibrium, one can capture the detailed fluctuation theorem even in an equilibrium 
simulation. This result is in contrast to many other fluctuation theorems which are valid 
far from equilibrium and at equilibrium they boil down to trivial identities as there is no 
average entropy production or consumption at equilibrium. 

Another issue is the choice of the microcanonical ensemble for the extended system. 
This is in contrast with Jarzynski's suggestion, which is to sample the reservoir degree of 
freedom from a Gaussian distribution. Instead, the present derivation treats all the degrees 



of freedom on equal footing and hence is more appealing, whereat 
are the usual ones of ergodicity and equal apriori probabilities 



re only assumptions made 

a. 



V. ILLUSTRATION OF DETAILED FLUCTUATION THEOREM IN DIFFERENT 
ENSEMBLES 

For a better appreciation of our results, we shall study the detailed fluctuation theorem 
in the context of a few popular extended phase space methods. All the examples mentioned 
below have the non-Hamiltonian evolution such that the system of interest evolves consis- 
tent with the appropriate ensemble the equations are supposed to mimic. The statistical 
mechanical properties of these systems are very extensively studied in the earlier papers 



[QQ 



15]] and hence our interest here will be limited to examining them in the light of the 



applicability of the detailed fluctuation theorem. 
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A. NOSE HOOVER DYNAMICS 

In the Nose- Hoover thermostatting scheme [3|, an N particle system in d spatial dimen- 
sions with Cartesian positions r = {ri,...,r]y} and momenta p = {p 15 ...,pn} interacting 
through N particle potential $(ri, ..., rjv) is augmented with 2 heat bath variables Y =(ri,p r] ) 
such that the 2dN + 2 dimensional extended phase space vector T = (r, p, rj,p v ) evolves as 
follows: 



n 



Pi 



Pi, 



III 



<9$({r;}) Prj 



dv 


i 


Q*" 


f) = 


Pn 

Q 




r n t,2 




- 




— 


iNk B T 


_i=i 







(35) 

(36) 

(37) 
(38) 



The parameter Q represents the strength of the Nose-Hoover coupling which controls the 
time scale over which the equilibration takes place and T is the temperature at which we 
wish to maintain the system of interest. A curious point that can be noted in these equations 
is that the variable rj does not explicitly get connected to other degrees of freedom. Still, it 
is profitable to retain this variable in the dynamical equations of motion as it facilitates the 
casting of these equations in the desired form (ITT]) 



V 

Pi 

\PvJ 



for 



A? 



10 

1 

-10 -pi 

-1 pi 



/ a$({rj}) \ 

dNk B T 



Ei 

m 

Pn 



n ^ = J2£r. + $ («) + %k + dNk * T ^ 



t=i 



'2m 



2Q 



Identifying the first two terms of the above as the system Hamiltonian, 



(39) 



(40) 



N 



^ S (W,{p}) = ^^ + $({r}) 



i=l 



(41) 
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we have 

2 

H{V) = H«({r}, {p}) + ^ + dNk B T V (42) 

Nose-Hoover dynamics is time reversal invariant and possess a "psuedo-Hamiltonian" as a 
constant of motion and hence, satisfies the essential requirements needed for the present 
analysis. The phase space compressibility of this system is given by eq ([1] 



2=1 i 

so the total phase space compression during an evolution from t=0 to t=r is from eq (fT5|) 

r 

AS = k B fA(T s )ds = k B dN(r](r) - r?(0)) (44) 

o 
There are several points in order. Given a time evolution of T = (r?, p^,f]°,p^) to T 4 = 
(r*, p;, r/ 4 ,p*), there is a phase space compression of ksdNirf — rj°). Note that this quantity 
is path independent and depends only on the initial and final state of the position variable 
of the heat bath, and does not explicitly depend upon the functional form of the potential. 
With this definition, one can gather all the trajectories that start at (za,^,P^) and evolve 
to (zb,?7 + AS,p v ) where ZAand zb are fixed initial and final states (positions and momenta) 
of the system of interest and i],p v a.nd p v are arbitrary (subject to the requirement that 
the extended state vector T stays on a constant energy hypersurface) . The probability of 
obtaining such a trajectory is P + (zb, -\-AS\za_). To obtain the anti-trajectories, from the 
time reversed final states one gathers all those states that start from (zg,^,^) and evolve 
into [z* A ,T] — AS,p v ). The probability of obtaining such a trajectory is P_(z^, — AS\z^). 

AS 

The detailed fluctuation theorem requires the ratio of these two probabilities go like e k s . 

Further, in the limit of Q — > oo, the equations of motion eqf )35p and eq( |36|) are reduced 
to that of a Hamiltonian evolution and we have, from eqf l37p . i](t) = r](0) and hence, from 
eq (T4"4"|) . AS" = 0. This means that the probabilities of forward and backward trajectories 
are identical: 

P(z A ~ z B ) = P(z* B ^ zl) (45) 

as can be expected for an isolated system. 
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Now, consider the probability that the system ip evolves from state ZAto state Zb in a 
time r. This probability, let it be denoted by P(z^ ►— > zb) , is 



P(z A ~ z B ) = / P(z B , A5|z A )rf(AS) (46) 

Now, applying the detailed fluctuation theorem on the right hand side, we have 



P(z A ►-» z B ) 



AN 



e^ P(z* A} -AS\z* B )d(AS) 



(47) 



For the Nose-Hoover thermostat, we know that AS = fee diV (77(7") — T)(Q)) and from eq 
we have 



7?(0) 



E-HHz A ) 



g|(0) 

2g 



^( r ) 



p 2 (r) 

P-^(ZB)-^ 



(48) 



diVfc B T ,v y dNk B T 

where E is the constant energy over which the microcanonical ensemble of the full system 
is defined. Hence we have, 



AS 



H s (z A ) 



PS(0) „., , , Pt(r) 



>B 



dN( V (r)-r](0)) 



2Q 



H s (zb) + 



1Q 



knT 



(49) 



Inserting this into the right hand side of eq (14"71) . 



h s (*a)- 



P(z A >-> z B ) 



2Q 



-ff s (z B )+" 



P?(t) 



1Q 



■p(z* A ,-AS\z* B )d(AS) 



(50) 



Since the z A and z B are independent of the integration variable, they can be pulled out of 
the integration and we have 



P(z A >-> z B ) 
Re-arranging the above, 



H S (z A )-H ( z Ei) 



P^(t)-p^(0) 



e 2 «^t P(z* A ,-AS\z B )d(AS) 



(51) 



-g S (^A) 



-g s (z B ) /" P%(t)-p%(0) 



(52) 



e fc s T P(z A ^z B ) = e fc s T /e 2 « fc s T P(z^, -AS|z*g)ci(AS) 

Assuming that the times t = and t = r are chosen sufficiently long time after the system 
has equilibrated, we can assume that the kinetic energy distribution is consistent with the 
equipartition theorem, and hence the total kinetic energy is equal to the ^dNksT. 
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N 2 1 

^t = 2 dNksT (53) 

With this assumption, and from the last of the Nose-Hoover equations eq( 155|) . we see that 
p v = 0, and hence e 2 Q k B T = 1. Substituting this in eq (1521) . we get 

e^^P(z A >-► z B ) = e^^ - / P(z^, -AS\z* B )d(AS) (54) 

The integral in right hand side of the above equation can readily be identified as the proba- 
bility P(z* B >—> z* A ) and since the system Hamiltonian also the time reversal invariance, this 
probability is also equal to P{zb >— > za). So, we have: 

e k B T P(z A >-> z B ) = e fc s T P(z B ^ z A ) (55) 

Which is nothing but the statement of detailed balance, which is valid for any system at 
equilibrium. 



B. NOSE HOOVER CHAIN DYNAMICS 

If more than one conservation laws are obeyed by the dynamical system, it is well known 
that the Nose-Hoover thermostatting scheme fails to generate the canonical ensemble. This 
happens because the accessible phase space gets restricted by the conservation and hence 
the system fails to access all the regions of the phase space in the course of its dynamics. 
This problem can be overcome by extending the number of heat bath variables. One such 
method is the Nose-Hoover chain thermostat. Its equations of motion are given by: 



m 



(56) 
Pi (57) 



<9r ; Qi 



% = ^p J = 1,-,M (58) 



Pm 



N 

"m 



V — - dNk B T 

' J rn ■ 



i=l 



P ^Pn, (59) 
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Pm 



P 



Vj-i 



Qi- 



-knT 



1\ 



0+1 



Qi 



3+1 



Pfb 3 



2,..,M-1 



(60) 



Pvm 



P 



Om-i 



knT 



Qm-i 
It can be shown that these equations of motion, preserve the psuedo-Hamiltonian: 



(61) 



N n 

pf 



N 1 

pl 



M 



n ^ = E^T + *(W) + E^T + WKTrn + k B T(J2v k ) 

i=l l 3=1 ^ fc=2 

The phase space compressibility for this system of equations is given by 



(62) 



A(r s ) = {dNfji + m + m + ... + m) 



(63) 



Again, the phase space compressibility rate for this system is also a total time derivative 
of the heat bath variables and the psuedo-Hamiltonian is invariant under time reversal. 
Hence we can derive the detailed fluctuation theorem in the context of Nose-Hoover chain 
thermostat. 

C. MTK ISOBARIC ENSEMBLE 



The Martyna, Tobias and Klein ensemble is also based on the Non-Hamiltonian phase 
space formulation and is known to generate the correct isobaric ensemble. The equations of 
motion read ( 15] ) 



Pi 



• Pi , Ps 

m W g 

0*({'i» p« 1 Tr(p g ) Vr 

"Pi Ar tjt Pi / ^> Pi 



<9r ; 



W n 



N, W g 



Qi 



Pg h 



p. =det[h](P^ - IP) + -L V ^Ll - ^p 



Qi 



(64) 
(65) 

(66) 
(67) 
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m = P f i* = f (68) 

Pvm = G m ( 69 ) 

Ph = Gj ~ ^Pfe P^m = Gm (70) 

where p( mi )is the internal pressure, / is the 3x3 identity matrix and G^are the thermostat 
forces, given by 

N n 2 P 2 

G 1 = y^- - dk B T G, = -^- k B T (71) 



1=1 



and 



6 X = Tr[ ^ Pg] - rf 2 A;T ^ = 4± - fc B T (72) 



These equations have the psuedo-Hamiltonian as 



* = E^+ $ ({ r })+^^+^( h )+E 



p% + p% 



2Qj 2Q 



M 



+N f k B T Vl +d 2 k B T£ 1 +k B Tj2(r ]k +tk) 



'2m,- VL J/ 2W, 

• j=2 L " J ~^JJ 

(73) 
The phase space compression rate for this system of equations comes out to be 

, M 

A = (d - l)-ln[det(h)] + dNfj! + d% + y \rj k + &] (74) 

fe=2 

It is readily evident that, as with Nose-Hoover, Nose-Hoover Chain and massive thermostat- 
ting schemes, the phase space compression rate is again a total time derivative of the heat 
bath variables alone and the psuedo-Hamiltonian is invariant under time reversal opera- 
tion. Further, the phase space compressibility is a function of position-like variable of the 
heat bath and hence satisfies the assumption, eq (j2"2"|) . This implies that the MTK isobaric 
ensemble is capable of capturing the detailed fluctuation theorem. 

Also, from eq (H3|) . eq ([63]) and eq( lT4"|) . we see that the phase space compression for 
all the extended system dynamics is dependent only on position-like variables of the heat 
bath and not on the system variables per se. This is evident from the extended phase space 
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fc=2 



formulation itself where there are no coupling terms between the heat bath variables and our 
Hamiltonian system variables. The phase space compression for the Hamiltonian systems is 
always zero, hence the contribution towards the phase space compression, eq [TH is from the 
heat bath variables alone. 

As is evident from the examples above, it is important to identify all the conservation 
laws satisfied by a given set of dynamical equations. For instance, in many thermostatting 
schemes, the psuedo-Hamiltonian, % is usually not the only conserved quantity. In such 
cases, the probability distribution has to be sampled from the hypersurface defined in eq 



a 



( TT9|) [3], We note here that all conserved quantities have to be time reversal invariant for the 
present proof to go through. 



D. GAUSSIAN ISOKINETIC ENSEMBLE 

Another example where the non-Hamiltonian phase space formalism can be readily ap- 
plied is the case of Gaussian isokinetic ensemble, which keeps the kinetic energy of the system 
constrained to a particular value but generates a canonical distribution in the coordinate 
space. The equations of motion read 



(75) 



Pi = Fi - 



Pi , 

•i = — I = 1 

mi 


,..,N 






" N 

EFj-Pj/mj 


Pi i = 


= 1,. 


.,N 


N 


EpIM 









(76) 



This isokinetic ensemble method different from other non-Hamiltonian phase space 
methods in that there are no extra degrees of freedom appended to the system, and also 
the total energy of the system is also not conserved. But, by construction, one has the 
conservation of the total kinetic energy and the unnormalized microcanonical probability 
density can still be defined as: 



TV 2 

p(p) = *(£^-#) 



j=l 



m 



(77) 
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where K is an arbitrary constant. The phase space compressibility rate of this system, can 
be obtained as 



N 



A = -^V ri .ri + V Pi .pi (78) 



i=l 

from eq (T75|) and eq (T76|) this becomes gives 

3iV-l#({r}) 

where the assumption is that there is no explicit time dependence of the potential on time, 
so that the partial derivative of the potential is zero. 
The Compressibility for this system is given as 

3JV-1 r fd6({r\) , 3N-l ri , , (/ .. 

A5 = —j -^-^ = -k B — — [0(r B ) - 0(r A )] (80) 

o 

The above equation implies, as the system starts from Za = (r A , Pa) and evolves to 

Zb = (rB?PB)) the entropy generated is proportional to the change in the potential energy 

at the end points, (J)(ta) — <P{ v b)- As the constraint, total kinetic energy, is invariant under 

time reversal and the phase space compressibility of the system eq flHOj) satisfies eq( | 22 |) . the 

present derivation of the detailed fluctuation theorem applies for the system evolving under 

the Gaussian isokinetic ensemble. 

The absence of external degrees of freedom in this example means that the phase 

space vector evolves deterministically and hence the probabilities P(zb, AS\za) and 

P(z B , — AS\z* A ) are reduced to just product of delta functions as there are no heat bath 

variables to integrate over. The joint probability P(zb, AS\za) is actually just a function of 

Zb and ZAalone as AS itself is a function of Zb and Za- In that case, consider the probability 

P(ta >— > i\b) = P{yb\ta) that the system evolves from a positionr^ to a position r^ in a 

time r. This involves integrating over all momentum variables and all possible phase space 

compression values. 



P(r A -» r B ) = f dp A dp B f P(AS)d(AS)P(z B = (r B , Pb), A5|z a = (r A , p A )) (81) 
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where P(AS) is the probability distribution of AS". Since the phase space compression 
depends only on the initial and final coordinates, eq (|80j) . we have 



P(AS) = S AS - AS{r A ,r B ) 



where 



AS(r A ,r B ) = -k B 



Inserting this into the above equation 



3N-1 f fd(f){r 



K 



dt 



dt 



(82) 



(83) 



P{v A ~ r B ) = / dp A dp B d(AS)5 AS - AS(r A , r B ) e^P(z* A , -AS\z* B ) (84) 
Now, AS'(r J 4, r B ) = —AS(r B , r A ) hence 



AS(r A ,r B ) 



dp A dp B / d(AS)S AS + AS{r B ,r A ) P(z* A ,-AS\z* B ) (85) 



P(r A >-> r B ) = e 



The RHS can be rearranged, remembering that the delta function is even in its arguments, 



P(r A >-> r B ) = e^ J d(-p A )d(-p B ) f d(-AS)S [AS - AS(r B , r A )] P(z% AS\z B ) 

(86) 
The integral on the right hand side can be readily identified as the probability of system 
evolving from r B to r^and hence 



Now, from eq(l83 



AS(r A ,r B ) 

P(r A >-» t b ) = e "b P(r B >-> r A ) 



AS{r A ,r B ) = -k B — — — [</>(r B ) - H t a)] 



(87) 



(88) 



Choosing the arbitrary constant K = (3iV — l)k B T, where T is the desired temperature, we 
have 



( r A)-<t>{r-B) 



P(v A ^v B ) = e k BT p( rs _» rA ) 



(89) 
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Rearranging the above equation, we see that the detailed fluctuation theorem just boils 
down to the detailed balance equation in the configuration space: 

0Oa) 0( r B) 

e "bt p(r A >-» r B ) = e fe s T P(r B >-» r A ) (90) 

The choice i^ = (3iV— l)ksT is natural from the equipartitioning theorem which says that 
every independent momentum degree of freedom has an average kinetic energy of \k%T . The 
system has (3iV — 1) independent momentum degrees of freedom (3iV momentum variables 
and one constraint on the total kinetic energy), hence the total kinetic energy is |(3iV — 
l)ksT. The Gaussian Isokinetic ensemble is a simple yet powerful example to appreciate that 
the Non-Hamiltonian characteristic alone is sufficient to bring out all the seemingly counter- 
intuitive features of the statistical mechanics like the entropy production and consumption 
anisotropy, direction of time emerging from the time reversible dynamics, anisotropy in the 
transition probabilities, emergence of Boltzmann distribution in the configuration space etc. 

VI. DETAILED FLUCTUATION THEOREM AND CONSERVATION LAWS 

Consider the Nose Hoover thermostatting scheme. If the forces acting on the system are 
derivable from a two body potential, 

N 

^(n,...,r JV ) = i^0 a (|r i -r i |) (91) 

such that the net force acting on the system is zero, then there are d additional conserved 
quantities that emerge in the Nose-Hoover dynamics: 

Pe v = K (92) 

n 

where P = J^Pj is the total momentum of the system and K is an arbitrary vector in d 
dimensions. In our context, this would mean that we can no longer use eqs ( I2"0"j) but instead 
use p(T) = 5(H (r) — E)5(Pe r} — K) as the correct distribution function. The problems 
associated with the presence of additional conservation laws in the context of generating 
the dynamics appropriate to a desired ensemble is well studied [3|, where it is shown that 
the presence of hidden conservation laws in the dynamics will lead to an ensemble different 
from the required canonical ensemble. It should be noted that the suggested solution of 
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appending an extended chain of thermostats, instead of one, though known to generate the 
correct distribution in the system subspace, will not invalidate the presence of additional 
conservation laws. So, the augmentation of the Nose-Hoover to Nose-Hoover chain will not 
help in demonstration of the detailed fluctuation theorem. The problem arises from the fact 
Pe v is not invariant under time reversal: (Pe v )* ^ Pe v and it is because of this problem, the 
present method of derivation hits a roadblock. In fact, for all systems which have conserved 
quantities which do not have a definite parity or are of odd parity under time reversal would 
fail to capture fluctuation theorems of the usual type, as the time reversed states are not 
accessible to the system. But for systems which have conserved quantities which are of odd 
parity under time reversal, sayK, it is easy to see that the Detailed Fluctuation Theorem 
takes the form: 

P + (z B} AS\z A} K) = e^P4z* B ,-AS\z* Al K*) (93) 

VII. FREE ENERGY RELATIONS FROM THE DETAILED FLUCTUATION 
THEOREM 

It would be worthwhile to investigate whether the free energy equality, Jarzynski's iden- 
tity, is derivable from the detailed fluctuation theorem in the non-Hamiltonian framework. 
For simplicity, we shall attempt to derive the Jarzynski's identity from the detailed fluc- 
tuation theorem result for a system coupled to Nose-Hoover thermostat. The Jarzynski's 
identity [18] reads, 

< exp(-/3W) >= exp(-(3AF AB ), = -L (94) 

where AF^is the equilibrium free energy difference between A and B: 

AF AB = -k B Tln f^-j (95) 

where Z A and Z B are canonical partition functions of the systems A and B. 

Z A = f e-^ H ^ rA ' PA) dr A dp A (96) 

and similarly 



21 



Z B = f e-P H ^ r *> p ^dr B dp B (97) 

where H^and H B are the Hamiltonian of the two systems whose free energy difference is 
to be computed. The Hamiltonians are superscripted to make the distinction between the 
Hamiltonian of the system of interest and the psuedo-Hamiltonian, which contains reservoir 
degrees of freedom also. We assume that there is a single time-dependent Hamiltonian, 
which at the time t = is the Hamiltonian corresponding to the state A, H A and at a 
time t = t transforms to the Hamiltonian corresponding to the state B, Hb- This variation 
can be brought about by, for example, a time dependent potential $(ri, ..., r/v,£), such at 
$(ri,...,rjv,0) = $A(ri, ...,rjv) and $(ri, ...,r N ,r) = $b(ti, ~.,r N ) where $ A and $# are 
the potentials of states A and B respectively. 

The change in the energy of the "system of interest" due to this time variation of the 

n n 

potential is given by [91. Illl|: 

T T 

H B -H A = H S ( ZB , t) - H 8 (z A , 0) = fdt dHa ^ t] z t + Jdt dHS ^ t] =Q + W (98) 

o o 

The above equation can be called the mathematical formulation of the First law of Thermo- 
dynamics, where the term on the left hand side is identified with the change in the internal 
energy of the system, the first term of the right is identified as the heat Q transferred from 
the bath to the system and the second term is the work performed on the the system. Note 
that both, work and heat are defined in terms of the system Hamiltonian alone. The effect 
of thermostatting is felt only through the coupling of evolution of the system variables and 
heat bath variables. 

From the explicit functional form for the system Hamiltonian and the Nose-Hoover equa- 
tions of motion, we can calculate the first term on the right hand side of eq (198]) as 



t_1 

V (t) N 2/ 

Identifying the term -y- — as i] and J^ Pi J^ - as p v + dNkgT from the Nose-Hoover equations 

Q\ i=l m 

of motion eq (38) and eq (39) we get 
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T 



Q = - fr](p v + dNk B T)dt (100) 



o 



As discussed earlier, If the times t — and t = r are such long after the system has reached 
steady state such that the average kinetic energy is determined by the temperature, we have, 
from eq (l38|) and eq fl53|) 



Q(t) = -dNk B T fdtr){t) = -dNk B T(r](r) - 77(0)) (101) 

'0 
Identifying dNi)(t) as the phase space compression rate for this system, eq( H3|) . we have 

T 

Q( T ) = k B T I \{t)dt = -TAS (102) 

'0 

Thus we see that k B T times the total phase space compression can be identified with the 

heat lost by the system to the thermostat. This is consistent with the assumption of eq Q). 

If the term AS" can be identified by the chan ge i n the entropy, then the above equation boils 



to the Second Law of Thermodynamics, 



16, 



ll| 



With this identification, we are ready to take-on the Jarzynski's work theorem in the 
context of Nose-Hoover thermostatted system. 



< eHWJ >: 



f dz A p(z A )P + (z B , AS\z A )e^ w ^ z »' As ^d(z B )d(AS) (103) 

It is easy to understand how the above integral is constructed, the quantity P + (z B , AS\z A ) 
gives the probability of system making a transition from z A at time t = to z B at time 
t = t and in the process generates a phase space compression of AS*. p{z A ) denotes the 
probability that the system is found in the state z^at time t — 0. If the times are so chosen 
that the system is fully equilibrated at time t=0, this probability is actually the probability 
of the canonical ensemble: 

e -/Jtfi(Z A ) 

PM = -^ (104) 

where 7i A is the partition function for the system in the state A, eq (I96p . Substituting this 
and also the detailed fluctuation theorem, we have 
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< e-? w >-- 



J- f dz A e- 0H ^ ZA) e^i P(z%-AS\z B )e^ w( - ZA ^ As ^d(z B )d(AS) (105) 

Li A J 



z 

Consider the terms in the exponentials, from eq (ll05l) and eq (198 j) we have 



and thus 



- [H s A {z A ) - TAS + W(z A , z B , AS)} = ~P{H s A {z B ) (106) 

< e ~ m >= ^- f dz A e-P H z {zB) P(z A , -AS\z* B )d(z B )d(AS) (107) 
Za J 

P(z* A , -AS\z B )dz A d(AS) = 1 (108) 



But 



And hence we have 

< e -PW >= _L f d ZBe -/»^(*B) = |£ = expt-pAFAB) (109) 

Za J Za 

This is the Jarzynski's Work theorem, eq fl95l) we set out to prove. 

Although the connection between the detailed Fluctuation theorem and the Jarzynski's 
identity has been established here only for the case of Nose-Hoover thermostatting scheme, 
it should be evident that the Jarzynski's work theorem can be derived in all contexts where 
the Detailed Fluctuation theorem is applicable. 

For illustration, consider the case of Gaussian Isokinetic ensemble. As mentioned above, 
these equations of motion fail to generate the proper canonical sampling in the momentum 
space but generates a canonical distribution in the coordinate space. But from eq, we see 
that the Free energy differences depends on the logarithm of the ratio of the two partition 
functions and hence the momentum partition function cancels out in the ratio and we are 
left with the ratio of the configuration partition functions at the states A and B. So we can 
see that one can realize the Jarzynski's identity in the Gaussian Isokinetic ensemble. 

As with the Nose-Hoover thermostat example above, Consider a system given by the 

N p 2 

Hamiltonian, H s ({r}, {p},£) = ^7r J - + ( ^ ) ({ r },^) we have the Gaussian Isokinetic equations 

i =1 2rrii 

of motion of the form, 

fi=— i = l....N (110) 

mi 
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Pi = -V r .0({r}, t) + a({r}, {p}, t)p, i = 1, .., N 



111) 



where 



«({r},{p},t) 



' AT 

j=l 

N 

Epf/ m j 

i =1 



(112) 



and the potential </>({r},£) is such that $(r l5 ..., rjy, 0) = $^(r l5 ..., rjy) and 
$(r l7 ..., r^v, t) = $s( r i) •••? r Ar) where $^and $# are the potentials of states A and B respec- 
tively. The time variation of this potential indicates that the work is done on the system. 
Consider again, eqf )98|) . we have the heat lost to the thermostat given by 



Q= dt 



d H s (z t ,t) 
dz t 





substituting for a({r}, {p},£) and simplifying we have 



f N 

z t = dta({r},{p},t)^2pi/ncii 



Q= dt 
o 



N 



^V rj 0({r},t). Pj -/mj 

Li=i 



AT r 

Identifying £ V rj 0({r}, t^/mj = ^A{t) we have Q = ^fdtA(t) 

3=1 



*J ~~ SAr- 
Together with eq f)44p this gives 



Q 



-K AS 



(113) 



(114) 



(115) 



3N -I k B 
Choosing the arbitrary constant K = (3N — l)fc#T (for reasons already discussed ) we have 



Q = -TAS 



(116) 



Again, we see that the heat lost by the system is proportional to the total phase space 
compression, as with the Nose-Hoover thermostat. With this identification, the procedure 
to calculate < e~^ w > is essentially unchanged from the Nose- Hoover case, and we have the 
Jarzynski's identity in the case of a system coupled to a Gaussian Isokinetic ensemble. 
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VIII. CONCLUSION 

Detailed Fluctuation Theorem has been extended to a class of thermostatted systems, 
evolving under the extended system dynamics. It is demonstrated that this theorem retains 
the same form as for the original DFT for entropy production when one replaces the thermo- 
dynamic entropy with phase space compressibility. This theorem is of a wider applicability 
than its original counterpart and can be applied even to the systems at equilibrium. It is 
shown that this detailed fluctuation theorem is formally equivalent to the detailed balance 
equation for systems at equilibrium. Rederivation of the Jarzynski's identity through the 
Detailed Fluctuation theorem has been demonstrated for both Nose- Hoover thermostat and 
the Gaussian Isokinetic ensembles. 
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